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Abstract 



We study a model of stochastic deposition-evaporation with re- 
combination, of three species of dimers on a line. This model is a gen- 
eralization of the model recently introduced by Barma et. al. (1993 
Phys. Rev. Lett. 70 1033) to q > 3 states per site. It has an infinite 
number of constants of motion, in addition to the infinity of conserva- 
tion laws of the original model which are encoded as the conservation 
of the irreducible string. We determine the number of dynamically 
disconnected sectors and their sizes in this model exactly. Using the 
additional symmetry we construct a class of exact eigenvectors of the 
stochastic matrix. The autocorrelation function decays with different 
powers of t in different sectors. We find that the spatial correlation 
function has an algebraic decay with exponent 3/2, in the sector cor- 
responding to the initial state in which all sites are in the same state. 
The dynamical exponent is nontrivial in this sector, and we estimate 
it numerically by exact diagonalization of the stochastic matrix for 
small sizes. We find that in this case z = 2.39 ± 0.05. 

PACS Numbers: 02.50+s, 75.10J, 82.20.M, 05.50+q 
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1 Introduction 



Recently a very interesting stochastic model with deposition and evaporation 
processes has been introduced by Barma et. al. |l|, |2|. In this model one 
deposits atoms on to k adjacent vacant sites of a (^-dimensional lattice and 
evaporates atoms from any k adjacent occupied sites, with specific rates for 
deposition and evaporation. The cases k = 1 and k = 2 are exactly solvable, 
the former being trivial, and the latter being equivalent to ferromagnetic 
Heisenberg spin-| chain. For k > 3, on a linear chain of length L, the phase 
space of this fc-mer model consisting of the 2 L possible configurations is found 
to break up into an exponentially large number of dynamically disconnected 
sectors [[[], 0. This may be understood as being due to the existence of 
an infinite number of independent conserved quantities in this model [[|]. 
These conservation laws also give rise to a wide range of relaxation behavior 
[U]. For example the density- density autocorrelation function shows different 
power law decays in different sectors (t _1//4 , t -1 ^ 2 , t~ ' 59 ) and even a stretched 
exponential decay in some sectors. 

For k > 3, though there exist an infinite number of constants of motion, 
these still have not enabled us to get a full solution of the model so far. In 
this case the number of sectors and their sizes in the steady state has been 
calculated exactly, but the dynamics is understood only qualitatively. The 
quantum Hamiltonian corresponding to the stochastic matrix has fc-body in- 
teraction terms which makes the problem difficult to tackle analytically. Thus 
it is worthwhile studying models which show the same qualitative features 
and are at the same time more tractable analytically. 

With this motivation we define a variant of the k-mer deposition-evaporation 
model with q states per site. This model for q = 2 corresponds to the dimer 
deposition evaporation model studied by Barma et.al. For q > 2 this model 
shares many qualitative features with the trimer model, like the existence 
of an infinity of conservation laws encoded as the conservation law of the 
"irreducible string" ||, and different power law decay of the autocorrelation 
function in different sectors. However our model has the two advantages: 
it has only 2-body interaction and it has an additional discrete symmetry 
group of very large order. These properties make this model more tractable 
analytically than the k-mer model. In fact, one can get many exact eigenvec- 
tors of the stochastic matrix in our model very easily. However, we have not 
succeeded so far in solving this model completely. For simplicity we study 
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only q = 3 case in this paper, for higher q results will be qualitatively the 
same. 

This paper is organized as follows. In section 2 we define our model 
and write the stochastic evolution operator in the form of the Hamiltonian 
of a quantum many-body system of 3 species of particles. In section 3 we 
show that the vector space consisting of 3 L configurations breaks up into 
an exponentially large number of dynamically disconnected sectors, and also 
calculate their sizes. In section 4 we describe the additional symmetries 
of this model and study the resulting further decomposition of the phase 
space formed by all the configurations belonging to a sector into subspaces. 
We construct a class of exact eigenvectors in the next section. In section 
6 we study the decay of autocorrelation function in various sectors. As in 
the case of trimer model ||, we find sector dependent decay. The spatial 
correlation function in the steady state for the sector corresponding to an 
initial condition of all sites in the same state is calculated in section 7. It 
has power law decay with exponent 3/2. A numerical diagonalization study 
of the stochastic matrix to find the dynamical exponent in this sector is 
presented in section 8. 

2 Definition of the model 

Our model is defined as follows: At each site of a one dimensional lattice, 
there is a spin variable qi which can be in any of the three states a, b or 
a These spin variables undergo a stochastic time evolution given by the 
following rule: any pair of adjacent spins, which are in the same state, can 
flip together to any of the other two states with some specified rates for the 
transitions. For example an aa pair can flip to become either a bb pair or 
a cc pair. This process can be thought of as the evaporation of a aa dimer 
and immediate deposition of a bb or cc dimer at that place. We shall call 
this model the dimer deposition-evaporation model (DDE model). This is a 
special case of the general reaction-diffusion process (a + b ^ c + d) recently 
studied by Dahmen 0. We consider only the case of equal rates for all the 
transitions in this paper. 

Any configuration C on the lattice can be represented by the basis vector 
\qi, q 2 , ■ ■ ■ , qi), in a vector space of dimension 3 L . If P(C, t) denote the prob- 
ability of the configuration C at time t, then the master equation describing 
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the evolution of these probabilities can be written as 

j t \P(t)) = W\P(t)) (1) 

where \P(t)) = ^ c P(C,t)\C) and W is the stochastic matrix. 

We can write W in the form of a quantum Hamiltonian as follows. Con- 
sider a quantum mechanical system of 3 species of particles on a one dimen- 
sional lattice of length L with periodic boundary condition. Each species is 
characterized by a color q which can be either a, b or c. Particles of the same 
colour experience a hard core interaction so that not more than one particle 
of the same color can occupy a given site. Therefor every site can be either 
empty or occupied by up to 3 particles of different colors. Thus there are 8 
possible states on every site. Let g« and qf denote the Pauli operators which 
annihilate and create a particle of color q at the site i. Now consider the 
Hamiltonian 

L 

H — ^2 A iji+1 + J3 i)i+1 + C^j+i (2) 
i=i 

where 

A iti+ i = (6j b i+1 + Cj c i+1 — 2aj a i+l )didi + i 
Bi,i+i = (atat +1 + cfci +l - 2btbt +1 )bib i+1 
= {atat+i + i>tH+i ~ 2c+CiVi)Qc i+ i 

where the operator A iti+1 acting on any state results a nonnull vector only 
when q,i = q i+ \ = a and is given by, 

. . , a, a, . . .) = |. . . , b, b, . . .) + |. . . , c, c, . . .) - 2\. . . , a, a, . . .) (3) 

The operators B i>i+ i and C iji+ i have similar actions. The operator 

hi = afcii + bfbi + cfdi (4) 

counts the number of particles at site i. Clearly 

[h t ,H} = 0, (5) 

and the number of particles is conserved at each site. If we restrict ourselves 
to the subspace of the full Hilbert space in which there is only one particle 
at every site, ie. hi = 1 for all i, corresponding to every configuration 
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\qi, q2, ■ ■ ■ , qi) m this subspace there is a unique configuration specified by 
the same L-string in the dimer problem. The action of H on \{qi\) is the 
same as that of W . Thus W is represented by the quantum Hamiltonian 
given by (0) acting on the subspace with hi = 1 for all i. 

Make the particle-hole transformation — > qf) on all odd numbered 
sites for all q. We define the current operator by 

Ju +1 = afa i+1 + bfb i+1 + cfc i+1 . (6) 

Then the Hamiltonian up to the addition of a constant can be written as 

H — ^2i,2i+1^2i,2i+l + ^2i,2i-1^2i,2i-l + ^ n<i,a n i+l,a- (7) 

The dimer model now correspond to H operating on a sector where hi = 1 
for all even numbered sites and hi = 2 for all odd numbered sites. The 
corresponding stochastic process is one in which there is a constant rate of 
exchange of any two particles at nearest neighbour sites. 

3 The Sector Decomposition of the Vector 
Space 

As the original trimer model, this model has an infinite number of indepen- 
dent constants of motion. These are described most simply in terms of the 
"irreducible string" (IS) @. We consider free boundary conditions for con- 
venience. Each configuration can be represented as a string of L characters, 
each character being one of a, b and c. From the string corresponding to a 
given configuration C, delete all adjacent pairs of aa, bb and cc. Each such 
deletion decreases the length of the string by 2. Repeat this operation on 
the resulting string until no more deletions are possible. This defines the IS 
corresponding to C. 

It is straightforward to extend the arguments of to our model, and 
show that the IS is a constant of motion. In addition, any two configurations 
having the same IS can be reached from each other. Thus the IS can be used 
to label the different sectors into which the vector space breaks up. 

Consider the set of all configurations in which no adjacent pair of sites 
are in the same state. These configurations will not evolve under the given 
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dynamics and are said to be fully jammed. For these configurations the length 
I of the irreducible string equals L. The total number of such configurations 
is easily seen to be equal to 3 x 2 L ~ 1 . Now consider sectors which are not 
totally jammed (/ < L). The number of sectors labelled by irreducible strings 
of length I will be the number of distinct irreducible strings of length I which 
is 3 x 2 l ~ 1 . The total number of sectors can be easily seen to be 2 i+1 — 1. 

Now we proceed to calculate the sizes of these sectors, i.e. the number of 
distinct configurations which belong to a particular sector. The size of a fully 
jammed sector is clearly 1. This is because a totally jammed configuration 
can not evolve and hence is the only member of that sector. Let D(IS, L) be 
the size of the sector labelled by irreducible string IS on a lattice of length 
L. We define the generating function 

oo 

G(IS,z)=J2D(IS,L)z L (8) 

L=l 

where the summation over L extends from I, the length of IS to infinity. 
Of special interest is the sector corresponding to IS of length zero. This we 
will call the null sector and denote its irreducible string by <fi. To compute 
G(4>, z), we introduce a formal series which is a sum of all strings of arbitrary 
length, which are made up of three letters a, b and c, and reducible to (j). 

G(<f>) = 4> + aa + bb + cc + aaaa + bbbb + cccc + abba + . . . (9) 

A string S is said to be decomposable, if the irreducible string IS(S) = <fi and 
it can be written as S1.S2 such that IS(Si) = IS(S2) = <fi- If Gi((f>) denote 
the sum of all indecomposable strings which are reducible to then G(4>) 
satisfies the equation 

e <*> = T^iw (10) 

Further if Gj denotes the the sum of all indecomposable strings which are 
reducible to and starting with the letter q then 

Gj((f>) =G a I + G b I + G c j (11) 

Gj is given by 

Gj = q . = 4 — • 1 (12) 

1 - (Gf + G b I + Gj- Gj) 
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To get the generating function we replace all occurences of a, 6, c by z in 
equation (|TJ). Then G% G\ and Gj all reduce to the same power series in z, 
call it g(z), given by 



z 2 



9{z) = T^iz) (13) 

This equation determine g(z) as an explicit function of z given by 



g(z) = (1 - Vl - 8z 2 )/4. (14) 
Using equations (|10D and (PD, we get 

G( ^ ) = T^b) (15) 

The growth of -D(0, L) for large L is determined by the singularities of g(z) 
nearest to the origin. This happens at z c = ±l/\/8. And since g has a square 
root singularity near z c , it is easily verified that 

D{<f>,L) ~8 L/2 L~ 3/2 (16) 

for large L. 

One can easily generalize this procedure to find the size of the sector 
when the irreducible string is of finite length /. Let IS = aia 2 where 
on = a, b or c. Then G(IS) (defined similarly as G(<f>)) is given by 

G(IS) = 1 a, I a 2 I -a,G(<P) (17) 

1 — H {aci) 1 — H (a.2) 1 — H (an) 

where H(a) = Gj + G\ + G\ — GJ. Replacing a, b, c by z as before we get 

G <'**> = [rrWi=W (18) 

where g(z) is given by equation (0) and I is the length of the irreducible 
string IS. Then the size of the sector characterised by the irreducible string 
IS is the coefficient of z L in (|T8f), and hence 



ww-hi*^- (19) 
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where integration is over a small circle encircling the origin. From this it 
is straightforward to find the asymptotic behavior of D(IS, L) by doing the 
integral using the saddle point approximation. Let A = lini£,_ ( . 00 Z/L. Then 

D(IS,L)~[k(\)] L , (20) 

where 

fc(A) = (1 - A)~ (1 ~ A)/2 (1 + A)-( 1+A V 2 2( 3 ~ A )/ 2 . (21) 

Thus, in general the size of a sector increases exponentially with L and the 
growth constant depends on the density of the irreducible string. From this 
one can recover the results for the null and totally jammed sectors by noting 
that fc(0) = and fc(l) = 1. 

4 Additional Symmetries 

In the previous section we discussed the decomposition of the 3 L dimensional 
configuration space of the model into 2 L+1 — 1 mutually disjoined sectors. 
This sector decomposition described in terms of the irreducible string IS is 
maximal in the sense that it is possible to go from any configuration to any 
other configuration in the same sector. However, this only helps us to block 
diagonalize W. In order to do better, we have to find additional operators 
that commute with W. These operators are not completely diagonal in the 
configuration basis. We now describe how to construct such operators using 
the additional symmetries of our model. These symmetries are not present 
in the original k-mer model. They are also not present in our model if all the 
transition rates are not equal. Using these symmetries we are able to reduce 
the task of diagonalizing the stochastic matrix with in a sector, by further 
block diagonalizing it. 

The deposition-evaporation model can be recast in the form of a gener- 
alized interface growth model, similar to the KPZ model, where the variable 
at each site is a 2 x 2 matrix instead of a scalar. This is seen as follows: 
we consider three non commuting matrices A(a),A(b) and A(c) such that 
A 2 (a) = A 2 (b) = A 2 (c) = Id, the identity matrix. A simple choice which 
satisfies these conditions is the unimodular matrices 

*>-(ii).^>-(;iiH-(;^) 

(22) 
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Alternately, we could choose them to be 577(2) unitary matrices. For a 
configuration \C) = \qi,q2, ■ ■ ■ ,Ql), we associate a 2 x 2 matrix h{C) to 
every site % given by 

I i (C)=A(q i )...A(q 2 )A(q 1 ). (23) 

It is easy to see that the stochastic evolution of the dimer model can be cast 
in terms of matrices J's by the following local evolution rule: if Jj_i = Ij+i, 
then Ii is reset randomly to any one of the three values A(q)Ii_i at a constant 
rate say 1. The conservation law of the IS in this representation becomes 
the simple statement that II is unchanged. 

A configuration of dimers is then equivalently characterized by specifying 
the matrix variables Ii at each site i. This matrix formulation brings out 
clearly the existence of a discrete group of infinite order in the model. If 
we work with free boundary conditions the set of allowed values of the of 
the matrices J's can be put in one to one correspondence with the sites of 
a 3-coordinated Bethe lattice. Identify one site of the Bethe lattice as the 
origin and associate the identity matrix to this site. We color the bonds of 
this Bethe lattice by three colors a, b, c such that any three bonds meeting 
at a site are all of different colors. The matrix I a corresponding to a site 
a of the Bethe lattice is given by the ordered product of A(q)s (q = a, b, c) 
along the unique path from the origin to a. A sequence {qi} , % = 1 to L 
specifies a configuration of the DDE model, and also a unique L-step path 
on the Bethe lattice which starts at the origin. As an example, the path 
corresponding to the configuration accacc is shown in figure |l|. Configurations 
in the null sector correspond to paths that return to the origin. Any pair 
of adjacent spins that are in the same state in the dimer configuration will 
corresponds to an immediate retraversal of a step on the Bethe lattice. The 
stochastic evolution in the DDE model leads to a stochastic evolution of this 
L-step polymer chain on the Bethe lattice. An elementary movement of this 
polymer chain dynamics is illustrated in figure [I], where the transition from 
the configuration acccacc to abbacc is shown. 

The operation of interchange of one branch of the Bethe lattice with 
another is a symmetry of W since this only corresponds to a recoloring of 
the bonds. This is illustrated in figure where we have considered only one 
branch starting from the origin. At every site a on the Bethe lattice define 
an operator P a which interchanges the two subtrees starting from that site 
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and going away from the origin. P a changes a configuration C = {ij} to 
C = {11} where 1[ — > f(Ii), where / is a matrix valued function of its 
argument. 

Clearly [W, P a ] = 0. We say that site (3 is descendant of site a if the 
unique path connecting (3 to the origin goes through a. From the figure 
it can be seen that [P a , Pp] = if neither a nor (3 is a descendant of the 
other. If (3 is a descendant of a then P a Pp = P 7 P a , where P 7 is the operator 
at site 7 which is obtained from site (3 by the action of P a . In figure |2] 
for example, [P 2 ,P 3 ] = 0. But [Pi,P 2 ] ^ since 2 is a descendant of 1. 
Instead they satisfy the relation P1P2 = P3P1. This symmetry comes from 
the geometrical symmetry of the underlying Bethe lattice. This in turn is 
related to the internal color symmetry in the model, and we will call it 
recoloring symmetry C. 

In addition to these one has a permutation symmetry of order 3 (S3) about 
the origin which corresponds to interchanging the 3 color labels globally. The 
set of operators {P a } together with the operators corresponding to 5*3 will 
constitute an infinite set of conserved quantities, which are however not all 
commuting with each other. 

One interesting consequence of the recoloring symmetry is that the spec- 
trum of W in a sector depends only on the length I of the irreducible string 
and not on its details. This can be seen as follows: For a general sector, the 
endpoints of all the paths on the Bethe lattice will be fixed. The shortest 
path connecting these two points corresponds to the irreducible string. All 
P Q s such that site a does not belong to this path leaves the irreducible string 
and hence the sector it labels invariant. On the other hand one can go from 
one irreducible string to any another one having the same length by a series 
of branch flip operations using the P a operators which are on the shortest 
path. Since each of these operations are symmetry operations of W, i.e. 
[P a , W] = for all a, the spectrum will be identical for all these sectors. 
This property is not present in the original /c-mer model. 

We can make use of these symmetries for further block diagonalisation of 
W with in a sector. This is achieved by breaking the state space spanned by 
all the configurations in a particular sector into subspaces. We illustrate this 
with a few examples. 

For convenience we will consider the case of null sector. First consider a 
lattice of length 2. There are only three configurations in this sector which 
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are \aa), \bb) and |cc). In this case we can get all the eigenstates of the 
Hamiltonian using 53. Consider the states 

= \aa) + \bb) + |cc) 
1^) = \aa) +u\bb) + u 2 \cc) (24) 
\ip 3 ) = \aa) + u 2 \bb) + u\cc). 

where u is the cube root of unity. These are the eigenstates of the operator 
corresponding to cyclic permutation of a, b and c (P a bc) with eigenvalues 1, to 
and uj 2 respectively. Since [P a bc, W] = these are also the eigenstates of 
W. The eigenvalues of W corresponding to these eigenvectors are 0,-3,-3. 
Thus in this case we were able to diagonalize W completely. But this is a 
rather trivial example. 

Next we consider the case of L = 4. There are 15 states in this sector. 
First we consider the symmetry operations about the sites 1, 2 and 3 as 
shown in figure |3|. Consider the states 

l^i) = 
1^2) = 
1^3) = 

These are the non-degenerate eigenstates with eigenvalue —1 of Pi, P2 and 
P3 respectively. It can be verified very easily that these are also eigenstates 
of W. This gives us a decomposition of the 15 dimensional vector space into 
subspaces of dimensions 1,1,1 and 12. Using the cyclic permutation symme- 
try at the origin, we can get a further decomposition of the 12 dimensional 



-j^(\abba) — \acca)) 

^(\bccb) - \baab)) (25) 
-k=(\caac) — \cbbc)). 
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subspace. For this consider the eigenstates of P a bc 
I0i) 

102) 

10s) 

104) 



Ixi) 

|X2) 
1X3) 
1X4) 



16) 
16) 
16) 
16) 

Since P a bc commutes with W, {\(fii)}, {\Xi)} an d {16)} form 3 invariant sub- 
spaces under the action of W. Thus we have got the following decomposition 
of the 15 dimensional vector space spanned by all the configurations in the 
null sector in to subspaces of dimensions. 

1 + 1 + 1 + 4 + 4 + 4 = 15 (27) 

We have carried out the same procedure for the case of L — 6. In this case 
the size of the null sector is 87. Using the P a operators, this can be broken 
into 18 subspaces out of which 12 are of size 1, 3 are of size 7 and another 3 
are of size 18. 

A further decomposition of these subspaces can be achieved using the 
symmetry of dynamics under reflection. Using both the color and reflection 
symmetries, the resulting form of the decomposition of the vector space for 
the case L = 6 is 

87 = 12 x1 + 3x3 + 3x4 + 3x6 + 3x 12. (28) 



= -^(\aaaa) + \bbbb) + \cccc)) 
= -^{\aabb) + \bbcc) + \ccaa)) 
= -^(\bbaa) + \ccbb) + \aacc)) 
= -^(\abba) + \acca) + \bccb) 
+ \baab) + \caac) + \cbbc)) 

= ^ ( | aaaa) + uj \ bbbb) + oj 2 \ cccc) ) 
= ^ ( | aabb) + uj \ bbcc) + u 2 \ ccaa) ) 
= -j^{\bbaa) + uo\ccbb) + u 2 \aacc)) (26) 
= ^g((|a66a) + \acca)) + u(\bccb) 
+ \baab)) +uj 2 (\caac) + \cbbc))) 

= ^ ( | aaaa) + uj 2 \ bbbb) + uj \ cccc) ) 
= ^ ( | aabb) + uj 2 \ bbcc) + uj \ ccaa) ) 
= -^(\bbaa) + uj 2 \ ccbb) + uj\aacc)) 
= -l=((|a66a) + \acca)) + uj 2 (\bccb) 
+ \baab)) + uj(\caac) + \cbbc))) 
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where m x n means m subspaces of size n. 

In principle this procedure can be carried out for a lattice of any length. 
For a linear chain of length L, sites up to distance L/2 from the origin on 
the Bethe lattice are reachable. The number of sites that are one step away 
from the boundary and inside will be 3 x 2^ L_4 ^ 2 (L > 4). As we have seen 
for L — 4, the antisymmetric eigenstates of the P operators on each of these 
sites will be an eigenstate of W . Thus one can easily get a large number 
of eigenstates in this model. By constructing the antisymmetric eigenstates 
of the P operators on sites which are interior to this, one can get invariant 
subspaces as we have seen in the case of L = 6 (the 3 sectors of size 7). So 
far we have not been able to get a general result for the total number of 
subspaces and their sizes within a sector as a result of this symmetry. But 
one can use this symmetry for the numerical diagonalization of the stochastic 
matrix for small sizes, where one can restrict oneself to one subspace (say the 
symmetric subspace where all the eigenvalues of the P operators are 1). This 
we have done to estimate the dynamical exponent in the null sector, where 
we have seen that these additional symmetries lead to substantial decrease 
in the size of the matrices to be diagonalized. These results are discussed in 
section 8 

5 Construction of some exact eigenvectors of 
the stochastic matrix 

In the previous section we have constructed some of the eigenvectors of the 
stochastic matrix in the null sector using the recoloring symmetry of the 
dynamics of L-step chain on the Bethe lattice. In this section we show that 
the same procedure can be used to construct eigenvectors in sectors which are 
almost totally jammed, i.e. the length of the irreducible string I = L — 2n, 
where n is a small integer. 

As we explained in the earlier section, the spectrum of W is the same 
for all the sectors whose irreducible strings are of same length. So we need 
to consider only one representative sector which we choose to be the one in 
which the irreducible string is abab.... Consider the case when I = L — 
2. We use periodic boundary condition for convenience. There are two 
types of configurations in this sector. In the first type there will be three 
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adjacent sites, say i,i + 1 and % + 2 which are of the same color (either aaa 
or 666 e.g. \. . . abbbab . . .)). Let these states be denoted by \4>{i)). In the 
second type there will be two adjacent sites, say i + 1 and i + 2 which are 
having color c. (e.g. |. . .abccab. . .). Let this type of states be denoted by 
|x(i)). In both these types site % + 1 and site i + 2 have the same color and 
thus constitute a reducible part. It can be very easily seen that under the 
deposition-evaporation dynamics, the position of the reducible part makes 
a nearest neighbour random walk. Hence the dynamics in this sector is 
equivalent to a random walk problem with two states per site {\4>{i)) and 
|x(i))). The eigenf unctions \ipk) and eigenvalues of the stochastic matrix 
can be found by solving this random walk problem. 

Let 4> k (i) and Xk(i) be defined by <p k (i) = {(f>(i)\ipk) and Xk(i) = (xWlV'fc)- 
They satisfy the following eigenvalue equations 

h<Pk(i) = 4>k{i ~ 1) + <f>k(i + 1) +Xk(i ~ 1) +Xk(i) ~ ^<f>k(i) 

(29) 

AfeXfe(i) = <t>k{i) + <t>k{i + 1) - 2xfe(i) 
The solution for this eigenvalue equations are given by 

\ k = (3 - cos(A;)) ± ^3 + cos 2 (k) (30) 

and 

Mn) = e* kn 

XM = ^(^(l + e*)) ^ 

where k = 2irm/L, m = 0, 1, . . . , L — 1. 

In the case of free boundary conditions some of the localized eigenvec- 
tors can be constructed directly with out solving this random walk problem. 
Consider the situation where the reducible part is near one of the boundaries. 
i.e. the states | ab . . . abaa) , | ab . . . abbb) and | ab . . . abcc) . It is easy to see that 
the state obtained by antisymmetrizing, \ip) = \ab. . . abaa) — \ab . . . abcc), is 
an eigenstate with eigenvalue — 1. One more eigenstate can be constructed 
by antisymmetrizing about the other end. 

Now consider the case of I — L — 4. In this case the reducible part is of 
length 4. A typical configuration in this sector will be |. . . abaaab . . . abf3(3ab . . .). 
The reducible parts aa and (3(3 perform random walk under the dynamics. 
When these random walkers come closer to each other they can form state like 
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|. . . abaaaaab . . .) or an intermediate state like |. . . abc/3/3cab . . .). There are 
two such intermediate states which corresponds to the two values (3 — a,b. 
The two random walkers can go into such intermediate states and remain 
there for a while by making transition between them. This can be consid- 
ered as a short range attractive interaction between the two walkers. So the 
dynamics in this sector is equivalent to a problem of two random walkers on 
a line having a short range attractive interaction. 

We can use the same procedure of antisymmetrization to get some more 
the eigenvectors. One type of eigenvector is obtained by antisymmetrizing 
about both ends simultaneously. 

= [\bbab . . . abaa) — \bbab . . . abcc)} — [\ccab . . . abaa) — \ccab . . . abcc)} 

(32) 

The second type of eigenvector is obtained by antisymmetrizing two in- 
termediate states |. . . abcaacab . . .) and |. . . abcbbcab . . .). 

|^2) — I- • ■ abcaacab . . .) — |. . . abcbbcab . . .) (33) 

Since there are L positions on the string, where abcf3f3cab can occur, for the 
case of periodic boundary conditions, there are L such eigenvectors. 

This can be easily generalized for the case where there are 2n random 
walkers (2n < L). A typical intermediate state formed by pair wise union of 
random walkers will have n number of substrings of the form . . . abcaacab .... 
The eigenvector is constructed by antisymmetrization of all these. The cor- 
responding eigenvalue will be — 3n. Since these n substrings can be arranged 
in roughly L C n ways, we can get so many eigenvectors. As an example, in 
the case of n = 2 one such eigenvector is given by 

= \. . . abcaacab . . . abcaacab . . .) — |. . . abcbbcab . . . abcaacab . . .) 
+\. . . abcbbcab . . . abcbbcab . . .) — |. . . abcaacab . . . abcbbcab . . .) 

(34) 

6 Autocorrelation functions 

For the trimer model, the existence of the infinite number of conservation 
laws lead to non-exponential relaxation in equilibrium M. It is found that 
the density-density autocorrelation function decay as a power law for large 
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times and the exponent is different in different sectors. In most sectors the 
decay is t _1//4 but is of the form t~ 1 / 2 in most of the sectors with periodic 
irreducible strings. In some special sectors with periodic irreducible string, 
stretched exponential behavior (exp(-y/i)) is found. In the null sector the 
autocorrelation function is found to decay as t~ 3 ^ 2z , where z « 2.5. This 
diversity of relaxational behaviour has been explained in terms of the hard 
core random walkers with conserved spins (HCRWCS) model ||, ||. It has 
been argued on qualitative grounds that these two models are in the same 
universality class. However, a strict proof of this proposition is not yet avail- 
able. 

An important assumption of the HCRWCS model is that the spin carried 
by each walker is conserved, but does not affect the diffusive motion of the 
walkers in any way. This is not strictly true in the trimer model. However, 
for our model, this property can be exactly established as a consequence 
of the recoloring symmetry. Thus the argument which predict the different 
behaviour of spin-spin autocorrelation function in different sectors is much 
cleaner for this model than for the trimer model studied earlier. We proceed 
to give the arguments which are an adaptation of the arguments for the trimer 
deposition evaporation model 0, |(|, in some detail. The time dependent spin- 
spin autocorrelation function in the steady state is defined by 

C^(t) = Probst) = (3, ft (0) = a) - Probfe(t) = /3)Prob(g,(0) = a). (35) 

where Prob(x) denote the probability of event x and qi(t) is the state of the 
spin at the zth site at time t. Each of a and (3 is either a, b or a Since a 
and (3 each takes 3 values, one can define 9 autocorrelation functions in this 
model. Among these, by time reversal symmetry C l a At) = Cl a (t) and also 
Y^pC^Jt) = 0. Thus there are only 3 independent correlation functions. 

As in the trimer model this autocorrelation functions shows sector de- 
pendent decay [Q. Let us consider the following five different representative 
sectors characterised by the irreducible strings 1). ababab . . ., 2). acbcacbc . . ., 
3). abcabc . . ., 4). random string and 5). 4> (Null sector), where . . . denotes 
repetition. The length of the irreducible string in the first four cases is taken 
to be L/2. 

We will first present a qualitative theory for the decay of autocorrela- 
tion functions and then compare its predictions with our Monte Carlo data. 
The behavior of the autocorrelation function can be understood in terms of 
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the random walk of the characters which constitute the irreducible string. 
Consider a sector labelled by an irreducible string aj, a.i ■ • • , on, with the 
length of the irreducible string I is a finite fraction of the total length L. 
Let X\,X2 ■ ■ ■ ,xi be the positions of the characters in a given configuration 
C which do not get deleted under the deletion algorithm of section 3. Then 
q Xi (t) — aj (i — 1 to /) at all times. We can think of xi,x 2 ■ ■ ■ ,Xi as the 
positions of I interacting random walkers with the nth random walker from 
the left carrying a color a n = q Xn with it. The positions of random walkers, 
{xi} will change in time as C changes. The walkers move either to the left or 
right but always remain on the same sublattice. Further, they do not cross 
each other. 

We have argued in section 4 that the spectrum of the stochastic matrix 
in this model is completely independent of the detail of the irreducible string 
string sequence {c^} and only depend up on the length of the sequence. 
However, the autocorrelation involves a weighted sum of correlations of 
which gives rise to different relaxational behavior in different sectors. 

The important contribution to the autocorrelation function comes from 
times when the site i is occupied by character from the irreducible string. We 
shall assume with out proof that the contribution coming from times when 
the site i is occupied by a reducible character is qualitatively similar. 

Since the hard core random walkers will remain on the same sublattice 
they were initially, we can perform a sublattice average of C l a n. Let Y denotes 
the sublattice which is either A (odd sites) or B (even sites). In HCRWCS 
approximation we can write the expression for the autocorrelation function 
as 

+oo 

Clp{t)= £ P(k,t)M r (k) (36) 

k=— oo 

where P(k,t) gives the joint probability of finding two particles at a given 
site, one at time and the other at time t such that the difference between 
their labels is 2k. M r (k) is a measure of the correlation of the colors in the 
irreducible string, and is given by 

M r (k) = J2S(q Xm ,a)6(q Xm+2k ,(3), (37) 

m 

where the summation is over all the labels of the particles which belong to 
the sublattice T. The probability P(k,t) can be obtained using the results 
of random walks of hard core particles ]7|, ||. |9| . 
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Note that P(k, t) depends only on the total length of the irreducible string 
and not on its details. While M r (k) depends only on the irreducible string 
and is independent of time. This separation of the sector dependence and 
time dependence is the crucial feature which allows calculation of C v a ^ in 
different cases. This separation is exact due to the recoloring symmetry in 
our model, but only approximate in the original trimer model. 

Equation (j36l) can be used to find the behaviour of autocorrelation func- 
tion. For example consider the sector labelled by the irreducible string 
abab .... On the A sublattice q Xm = a for all m and hence 

+00 

<?£(*)= E p (k,t) (38) 

k=— 00 

This is the density- density correlation function of hard core particles which 
decays as t™ 1 / 2 for large t J7j. Our simulation gives the same power law 
decay which is shown in figure [|. From the symmetry between b and c on 
this sublattice and the constraint J2pC^Jt) = 0, it can be easily seen that 
C£g(t) ~ t~ 1 / 2 for all values of a and (3. The behavior of the correlation func- 
tions on both the sublattices is the same which follows from the symmetry 
of the irreducible string. 

In the case of the sector labelled by the irreducible string acbc . . ., M (k) 
is of the form c\ + Ci{— l) fe , where c\ and c% are constants independent of 
k. Thus C^ a (t) has two parts. The first part decay like t -1 / 2 as explained 
earlier. The second part is the Fourier transform of P(k,t) which goes as 
exp (-y/i) for large t. The stretched exponential decay dominates the short 
term behavior but asymptotically the behavior will cross over to r 1 / 2 . By 
symmetry C^ b {t) should have the same decay. On the B sublattice q Xm is c 
for all m, so the decay will be purely diffusive (t~ 1 / 2 ) as explained earlier. 
These behaviours can be seen in figure |5|. 

For the sector abcabc. . ., the decay of C^ a {t) is shown if figure §. The 
leading behavior of the correlation function can be shown to be stretched 
exponential decay by a similar argument as in the previous case. By sym- 
metry of the irreducible string the behaviour should be the same for all the 
correlation functions and on both the sublattices. This has been confirmed 
in our simulation. 

When the irreducible string is random, M r (k) is significant only for small 
values of k. For small k, P(k,t) ~ t -1 / 4 for large t 0]. Therefore the 
autocorrelation has a t" 1 / 4 decay which can be seen in figure 0. 
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Figure |7| also shows the decay of C^ a (t) in the null sector. It is a power 
law decay with the value of the exponent ~ 0.59. The same power law decay 
found for other values of a and (3 and on both sublattices, which follows from 
symmetry. Barma et.al. have obtained same estimate for this exponent for 
the trimer model in the null sector |4]] . This is an evidence that the dynamics 
in these two models for the case of the initial condition in which all sites are 
in the same state, belongs to the same universality class. We can understand 
the behaviour of the correlation functions in sectors where the irreducible 
string has a finite density using the HCRWCS model, this approach is no 
longer useful when studying the null sector. An analytical understanding of 
the temporal decay of autocorrelation functions in this sector is still lacking. 



7 Spatial Correlation Function in the Steady 
State 

Let C a ^{r) be the probability that, in the steady state, spins at two sites 
which are separated by a distance r are a and (3. We calculate this equal 
time spatial correlation function for the null sector. Let us assume periodic 
boundary conditions. Let 5* be the string corresponding to a given config- 
uration on this lattice with sites i and i + r are occupied by a and (3. Let 
Si be the substring of length r — 1 formed by spins between a and (3 and 5*2 
be the substring of length L — r + 1 formed by the remaining spins . The 
endpoints of the substring S 2 are a and (3. Let JS(Si) and IS(S 2 ) be the 
irreducible strings corresponding to Si and S 2 . Since IS(S) = <fi, IS(S 2 ) will 
be IS(Si) reflected about one end. We may write IS(S 2 ) = IS(Si)^ 1 . 

° aAr) ~ D(^L) (39) 

where the summation is over all the distinct irreducible strings of Si. D(IS, r— 
1) is the size of the sector labeled by IS on a lattice of length r — 1 with 
free boundary conditions and D^^IS^ 1 , L — r + 1) is the size of the sector 
labeled by IS -1 on a lattice of length L — r + 1 with boundary sites fixed at 
a and (3. 

Though we can calculate exact expressions for each of the above quan- 
tities, we derive the asymptotic behaviour by taking L — >• oo,r << L. The 
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length of the irreducible string I < r, therefor in this limit one can make the 
approximation D a ^(IS, L — r + 1) rs D a> p((f>, L — r — I + 1). In the case of free 
boundary condition the size of a sector depends only on the length of the ir- 



reducible string (|20| , |21|) and we may replace the sum over distinct irreducible 



strings by a sum over the lengths of the irreducible strings as follows. 

r M _ d mm, r - i)DaM L- r -i + i) 

CaAr) ~ D(^L) (40) 

Where N(l) is the number of distinct irreducible strings of length I, which is 
given by 3 x 2' -1 . For large L, D a> p((j),L-r-l + l) ~ (2V2) L ~ r ~ l+1 / (L - r - 
/ + l) 3 / 2 . Hence in the summation over I only terms corresponding to small I 
will contribute. We can find the asymptotic r dependence by approximating 
the sum by only the / = term. Then we have 

C a ^{r) ~ D((/>, r - l)D aiP (<j>, L-r + 1)/D(<j>, L). (41) 

For large L gives 

f(a,P) 



r 



for large r, where f(at,/3) is a constant independent of r but depends on a 
and (3. From the symmetry between different colors, it follows that f(a, (3) = 
k (1 - 35a,p) 



8 Numerical Diagonalization of the Stochas- 
tic Matrix in the Null Sector. 

As mentioned in the section 6, the decay of autocorrelation function in the 
null sector suggests that the dynamics in this sector belongs to a new univer- 
sality class. The description of dynamics in terms of the random walk of the 
characters of the irreducible string is not possible in this sector. One way of 
studying this sector is to directly diagonalize the stochastic matrix for small 
sizes. However this task is not very easy because the size of the matrix to 
be diagonalized grows very fast with the size of the lattice (~ (2v / 2) i / L 3 / 2 ). 
One can make use of various symmetries of the model to block diagonalize 
the matrix first and hence reduce the size of the matrix to be diagonalized. 
By using translation, reflection and the recoloring symmetry described in 
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L 


N 


N c 


Nr,t,c 


A L 


zl 


2 


3 


1 


1 






4 


15 


3 


2 


-10.00000 




6 


87 


12 


4 


-3.108548 


2.881701 


8 


543 


55 


10 


-1.629235 


2.245691 


10 


3543 


271 


26 


-0.974243 


2.304369 


12 


23823 


1399 


93 


-0.632500 


2.369330 


14 


163719 


7470 


338 


-0.438272 


2.379750 


16 


1143999 


40931 


1474 


-0.318921 


2.380690 


18 


8099511 


228918 


6801 


-0.240934 


2.380807 


20 


57959535 


1301778 


33746 


-0.187477 


2.381040 



section 4 we have been able to diagonalize the stochastic matrix for lattice 
sizes up to 20, in the fully symmetric subspace. The algorithm used here for 
diagonalizing the matrix using the various symmetries is similar to the one 



we have used for studying trimer model |10| . For details of the algorithm we 
refer the reader to this paper. 

In table 1, we have listed the size of the matrix (Nn t T,c) obtained after 
using all the three symmetries and the value of the largest nonzero eigenvalue 
in the symmetric subspace in the null sector for lattice sizes L from 2 to 20. 
For comparison we also list the total size (N) of the null sector, and also 
the reduced matrix size (Nc) if only the recoloring symmetry C is taken into 
account. We see that the additional recoloring symmetry gives a significant 
reduction in the size of the matrix to be diagonalized. 

The largest eigenvalue of the stochastic matrix is 0. The difference be- 
tween the largest and the second largest eigenvalues gives the gap in the 
spectrum of the stochastic matrix. We have computed this gap for lattice 
sizes from 2 to 20. By using the finite size scaling relation Al ~ L~ z , the 
dynamical exponent zi for a lattice of size L is estimated from the relation 

^ = log(-^-)/log(^). (43) 

The estimate of zl for various values of L are shown in table 1. The con- 
vergent of zl seemed to be quite good, and the estimate of the extrapolated 
value of z corresponding to L = oo is 2.39 ± 0.05. The error bar reflects our 
subjective estimate of the possible systematic errors in the extrapolation. 
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This value is in good agreement with the value of z = 2.5 ± 0.15 for the 



trimer model in the null sector [jig] , and is evidence that both these models 
are in the same universality class. 

We thank M. Barma for introducing us to the model, and for his contin- 
ued interest and discussions throughout the course of this work. We have 
benefited from several discussions with Peter. B. Thomas in the earlier stages 
of this project. 
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Figure 1: Dimer deposition-evaporation dynamics as a dynamics of polymer 
chain on the Bethe lattice. The transition between the states \accacc) and 
\abbacc) is shown. 



Figure 2: The operation of interchange of two branches that start from the 
site 3. This is equivalent to recoloring of the bonds, and is a symmetry of 
W. 



Figure 3: Construction of exact eigenvectors using recoloring symmetry. For 
L — 4, three exact eigenvectors can be constructed by antisymmetrizing 
about the sites 1, 2 and 3. 



Figure 4: The autocorrelation function C^ a (t) for the sector abab ... (o). The 
dashed line is a plot of t~ l l 2 . 



Figure 5: The autocorrelation functions C^ a {t) (o) and C^ b {t) (+) for the 
sector acbc. . .. For small t, C^ a (t) has a stretched exponential decay of the 
form exp(-t -1 / 2 ) but at later times it crosses over to pure diffusive decay. On 
the other hand C^ a (t) has pure diffusive decay at all times. 



Figure 6: The autocorrelation function C^ a (t) for the sector abcabc. . .. The 
decay is stretched exponential of the form exp(-t -1 / 2 ) which shown by the 
dashed line. 



Figure 7: The decay of C^ a (t) (o)in the sector where irreducible string is 
random, and in the null sector (+). These decays are approximately given 
by t -1 / 4 (dashed line), and t~ - 59 (dotted line). 
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